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Heat has always been a killing matter for traditional semiconductor machines. The underlining physical 
reason is that the intrinsic carrier density of a device made from a traditional semiconductor material 
increases very fast with a rising temperature. Once reaching a temperature, the density surpasses the 
chemical doping or gating effect, any p-n junction or transistor made from the semiconductor will fail to 
function. Here, we measure the intrinsic Fermi level {\Ep\ = 2.93 ic^T) or intrinsic carrier density {iijjj = 3.87 
X 10^ cm~^K~^- V), carrier drift velocity, and G mode phonon energy of graphene devices and their 
temperature dependencies up to 2400 K. Our results show intrinsic carrier density of graphene is an order of 
magnitude less sensitive to temperature than those of Si or Ge, and reveal the great potentials of graphene as 
a material for high temperature devices. We also observe a linear decline of saturation drift velocity with 
increasing temperature, and identifj^ the temperature coefficients of the intrinsic G mode phonon energy. 
Above knowledge is vital in understanding the physical phenomena of graphene under high power or high 
temperature. 

Since the discoveries of carbon nanotube and graphene^'^, a significant amount of attention has been directed 
toward the room temperature and low temperature measurements in pursuing a high carrier mobility^"^. As 
people seek device applications for graphene, high field and high current measurements start to appear^"^°. 
The importance of electron-phonon coupling is well recognized^ ^"^l However, little attention has been made to 
the role of intrinsic carriers in graphene, because the great ability of graphene to tolerate the impact of temper- 
ature masks the existence of its intrinsic carriers. Thus, the importance of graphene as a material for high 
temperature devices is not noticed. Here we use a Raman method to measure the level of intrinsic carriers in a 
voltage biased graphene device up to 2400 K. We find the level of intrinsic carriers in graphene is an order of 
magnitude less sensitive to the temperature than that of traditional semiconductors, silicon and Germanium. The 
result tells us that a properly fabricated graphene p-n junction or transistor can operate well above 1500 K and 
will not fail until other parts melt. Working together with silicon carbide (SiC), graphene might be one of the 
building materials for coming high temperature devices. We also observe a linear decline of saturation drift 
velocity in graphene with increasing temperature in contrast to a constant saturation drift velocity in traditional 
semiconductors, as the result of the degeneracy of carriers and the Pauli exclusion principle. At the same time, we 
identify the temperature coefficients of the intrinsic G mode lattice vibration energy up to 2400 K. Above 
information highlights the great application potentials for graphene in high temperature devices, suggests a 
new exploring direction for graphene research, and is critical in the physics of graphene under high field, high 
current or high temperature conditions. 

The graphene devices were fabricated by mechanical exfoliation from graphite. Figure la shows a sketch of the 
layout of graphene devices and electrical measurement. Contacts were formed on top of the graphene after steps 
of E-beam lithography, metal evaporation and liftoff. The contacts are 100 nm Au with 5 nm Cr between Au and 
graphene. The substrate is doped Si with 300 nm Si02 isolation layer on top. When an electrical bias is applied to a 
device, the drain electrode and Si substrate are grounded. A Kethley 2400 is used both as a voltage source and a 
current meter. During the measurements, the devices were kept in a flowing nitrogen gas environment in order to 
prevent the oxidization of the sample. A Horiba Raman microscope was used to measure the Raman spectra of the 
sample from the top. The excitation laser is a He-Ne 633 nm (1.96 eV) laser. The focusing objective is a long 
working distance SOX objective lens. Two graphene devices are measured; the optical pictures of these graphene 
devices are showing in Fig. lb and Ic. Device A is 5 jim in length and 2.7 jim in width, and device B is 4.85 jim in 
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Figure 1 | Graphene device information, (a), The sketch of the stacking structure of the devices and the electrical measurement layout, (b), Optical 
picture of the graphene device A. (c), Optical picture of the graphene device B. The widths and lengths of the devices are marked in the pictures. The blue 
dots are the optical focus positions during the measurements of the devices. The letters D and S indicate the drain electrodes and source electrodes. 



length and 17.1 jim in width. For example, we measured the device A 
from V^^ = 0 V to V^^ = 21 V with an increasing step of 1 V, then a 
downward scan (with a step size of — 1 V) back to Vsd = 0. For each 
different Vsd, ^ Raman spectrum from the center of the device was 
recorded; the currents of the device before and after each Raman 
measurement were recorded as well. 

Our measurement result reveals the dominating impact of 
intrinsic carrier due to increasing temperature in the devices under 
high bias conditions. In Figure 2, we plot the measurement results of 
Raman spectra as a function of Vsd for device A. We can see two main 
Raman peaks, G mode around 1600 cm"^ and 2D mode^^ around 
2600 cm~\ In both subfigures, G mode intensity increases with a 
larger Vsd- At the same time, a dramatic increasing of the background 
signal appears with a high Vsd voltage and this background will drop 
back down with the decreasing Vsd during the downward scan. In an 
earlier report^^, C.-F. Chen et al. presented that the G mode intensity 
has resonant behavior when the twice of the Femi level in graphene 
matches the excitation photon energy or the photon energy minus a 
G mode phonon energy. They also reported an increasing back- 
ground in the spectra near the resonant condition, and they attrib- 
uted this high background signal to hot luminescence in graphene. 
Our data's resonant behaviors in the G mode intensity and back- 
ground are very similar to those reported behaviors, therefore we 
believe that the G mode intensity change is due to a changing Femi 
level in the graphene device. The Femi level (Ep) in graphene is 
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Figure 2 | Raman spectra of device A plotted in color maps for different 

(b), The spectra during the Vsd downward scan. 



Vp is the graphene Fermi velocity, H is the reduced Planck constant. A 
higher Fermi level will have a higher carrier density, and vice versa. 
There are four possible reasons may change the device carrier den- 
sity. They are intrinsic carrier density increasing due to rising tem- 
perature, permanent defects in graphene, back gating effect due to 
increasing Vsd, and the desorbing of the surface molecules. In com- 
parison with Fig. 3a in C.-F. Chen et al.'s paper^^ after shifting the 
resonant peak due to different excitation photon energies (785 nm to 
633 nm, 1.58 eV to 1.96 eV), we use the relative intensity of the G 
mode at V^^ = 21 V over base intensity at Vsd ^ 0-5 V, and deter- 
mine that \Ep\ ^ 0.8 eV or n ^ 4 X 10^' cm"^ when Vsd = 21 V. First 
of all, this high carrier density clearly is not the result of permanent 
defects in the device, because the Vsd downward scan shows this 
carrier density will go down with a lower voltage. For unintentional 
residual and surface adsorbate doping, the carrier density is normally 
at the scale of 10^^-10^^ cm"^, and this doping should only decrease 
with an increasing temperature. The back gating effect of the Vsd can 
be estimated after calculating the back gating capacitance. For a 
300 nm Si02 layer (permittivity 8 = 4), the typical back gating capa- 
citance Cijg = 1.2 X 10"^ F-cm"^; for Vsd — 20 V, back gating carrier 
density at the middle of the device Ci,gVsd/2 = 7.5 X 10^^ cm"^. 
This value is two orders of magnitude smaller than what we have (n 
~ 4 X 10^^ cm"^). The only reasonable explanation for this high 
rising carrier density is the increasing intrinsic carrier density due 
to rising temperature within the device. Also, the scales of other 
effects, back gating and desorption of surface adsorbate, are a couple 
orders of magnitude smaller than the intrinsic carrier density under 
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Figure 3 | Measurement and fitting results for G mode intensities as a function of the V^j. (a), Device A during V^d upward scan, (b), Device A 
during Vsd downward scan, (c), Device B during Vsd upward scan. The data points are fitted with the curve from Fig. 3a in C.-F. Chen et al.'s 
paper^^ after shifting 0.38 eV (the difference between excitation photon energies, 785 nm to 633 nm or 1.58 eV to 1.96 eV). We only use the Vsd ^ 9 V 
data points for the fitting in order to give a better fitting weight for the high bias data points. 



the high bias condition in our sample. We determine that the ther- 
mally generated intrinsic carrier is the dominating factor for those 
short samples under high bias voltage conditions. 

First, we will use the G mode intensity to determine the Fermi level 
in the devices. The Fermi level in the graphene device is changed with 
Vsd, but not directly tuned by Vsd- The Vsd drives the graphene device 
to generate heat; the heat generated carrier density rises in the device 
with an increasing temperature; with the help of external field from 
Vsd, heat generated carriers break balance between electron carriers 
and hole carriers, and the major carrier winning out raises the Fermi 
level in the graphene device. For graphene with Fermi level at the 
Dirac point and no external field, Tian Fang et al. calculated ther- 
mally generated carrier density^^, rith = - ( 7^ | , where ks is the 

6 \nvp J 

Boltzmann constant, T is the temperature. Their result gives out a Uth 
for both the electrons and the holes, the electrons and holes cancel 
each other so the Fermi level is still zero. Clearly, our graphene 
devices under high field don't fit into their calculation's pre-con- 
dition, and their calculation result probably is numerically invalid 
for our case as well. In a connected real device with a certain electric 
voltage level, the Fermi level can go well above zero. For our devices, 
we assume the Fermi level due to rising intrinsic carrier density at the 
measurement point of a device, 

\Ep\ = ^Vs4+Eo, (1) 

where a is a device and location dependent constant. Eg is a constant 
Fermi level when Vsd ^ 0- The reasons behind this assumption are: it 
is reasonable to guess \Ep\(xkBT, and we notice that device temper- 
ature Thas a linear relationship with Vsd. (This will be discussed later 
when we discuss transport data, and in supplementary information.) 
Equation (1) works much better for Vsd downward scan than Vsd 
upward scan. During Vsd upward scan, the graphene device goes 
through desorption of surface adsorbate and increasing of perman- 
ent defects due to rising temperature^^, so the Eg is not a constant. 
During Vsd downward scan, the surface adsorbate was fully released 
at the upward scan step, and the permanent defects should be the 
same because the device is cooling, so the Eg is about a constant 
during downward scan. This argument is supported by Figure 4a, 
which shows that the downward scan data for G mode positions 
shows a much better smoothness and data continuity than the 
upward scan data. In addition. Figure 4a shows that the G mode peak 
position, Qq = 1590.44 cm~\ for Vsd ^ 0 after the downward scan. I. 
Calizo et al. reported that the graphene G mode energy linear coef- 
ficient with temperature is —0.016 cm"VK for a temperature range 
83 K-373 K, and Qg = 1584 cm"^ for T = 0 so Qg = 



1579.2 cm"^ for T = 298.15 K (room temperature) if the Fermi level 
in graphene is zero. The blue shifted 1 1.24 cm~\ in G mode energy at 
Vsd — 0 after downward scan, is due to the Fermi level Eg in graphene. 
Because^' Ar2G = \Ep\ X 42 cm'^eV^ we willhaveEo = 0.2676 eV 
for device A during downward scan. Figure 3 shows our measured G 
mode intensities as a function of the Vsd- The fitting curves are the 
curve retrieved from C.-F. Chen et al.'s paper^^ Fig. 3a for Ig vs. 2\Ep\, 
after blue shifted 0.38 eV. (They use a 785 nm-1.58 eV excitation 
laser and we use a 633 nm-1.96 eV excitation laser.) Our experi- 
mental data are fitted to the curve with two fitting parameters, a in 
equation (1) with constant Eg = 0.2676 eV and a linear scaling factor 
between our G mode intensity and Ig in C.-F. Chen et al.'s paper. The 
fitting is most accurate for device A during downward scan, because 
the Eg value is measured for it. We do the same fitting for device A 
upward scan anyway, just to show a comparison. We only use the Vsd 
> 9 V data points for the fitting in order to give a better fitting weight 
for the high bias data points, which matter the most in this fitting. 
The G mode intensities at the high bias voltages already rise signifi- 
cantly above the random error range, especially for the data from 
device A. Therefore, our data have enough signal-to-noise ratios to 
make this fitting method valid, although our data doesn't cover the 
full resonant peak. For device A, the upward scan data fitting gives a 
= 0.0520 eV/V; the downward scan data fitting gives a = 0.0513 eV/ 
V; these two values are close and we take the downward scan a = 
0.0513 eV/V as the value for device A. Data from device B has lower 
quality than the data from device A, however we use device B for the 
verification purpose only. Because device B broke down after the 
measurement at Vsd = 17 V during the upward scan, no data is 
available for device B Vsd downward scan and we can't accurately 
determine the Eg. In Fig. 4a, the high bias {Vsd ^ 6 V) data points' 
behaviors are similar between the upward scan and downward scan 
for device A, this might suggest that the average Eg value for high bias 
conditions in upward scan is close to the Eg value in downward scan. 
We also consider device B has a similar length and goes through 
similar processes as device A, it is likely that device B has a similar 
defect level as device A. So, we use Eg = 0.2676 eV from device A, as 
an approximation value, for the fitting of device B high bias data in 
Fig. 3c. The fitting gives a ~ 0.0566 eV/V for device B. 

Next, we will find out the temperatures at the measurement point 
for difference Vsd situations using G mode peak position informa- 
tion. The key information we used to determine temperature is: the 
value of room temperature, the intrinsic G mode peak position shift 
(after correcting the doping effect) as a result of temperature change, 
the reported G mode energy temperature coefficient for temperature 
range 83 K to 373 K from Ref 18, and a 2""^ order polynomial fitting 
of our G mode peak position. The 2"^^ order fitting gives the relation- 
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Figure 4 | The G mode peak positions as a function of the V^j. (a), Device A. The blue square, red cross, and blue circle data points are the 
experimental data points for V^^ downward scan, V^d upward scan, and the G mode position for downward scan after taking out the effect of carrier 
density, respectively. The blue curve is the fitting result for the corrected G mode position with a 2"^^ order polynomial function, (b). Device B. The blue 
square and blue circle data points are the experimental data points for V^j upward scan and the G mode position after taking out the effect of carrier 
density. The blue solid curve and dash curve are the fitting results for the corrected G mode position and the original data points, respectively, with a 
2"^^ order polynomial function. Only the data points of V^^ > 6 V are used for the fittings of device B. 



ship between the first and second temperature coefficients; our 
device temperature range overlaps with Ref. 18, so the reported the 
linear temperature coefficient from Ref. 18 sets another condition for 
the first and second temperature coefficients. Combining both con- 
ditions, we can have the values of the first and second temperature 
coefficients for G mode energy. With the value of room temperature 
and G mode Raman shift of each data point, we can have the device 
temperature of each data point. Figure 4 shows the G mode peak 
positions as a function of the Vsd for device A and device B. In 
comparison between upward scan and downward scan G mode posi- 
tions for device A, the downward scan data shows a much better 
smoothness and data continuity. The surface adsorbate was fully 
released during upward scan, so no chemical doping changes during 
downward scan. The device goes through cooling during downward 
scan, the permanent defects have a constant amount, and only shift 
the downward scan data a constant amount in G mode energy. 
Therefore, we use the downward scan data for further data analyzing. 
During Vsd downward scan, the G mode peak position, Qq, satisfies 
Qg = ^G_T + AI2g /n + AQG_defect^ where Qg_t is the intrinsic G 
mode energy as a function of the temperature, AQG_in is the G mode 
energy shift due to rising Fermi level (intrinsic carrier density) as a 
function of Vsd, AQG_defect is a constant from the permanent defects 
level in the graphene. We know that AQgju + AQG_defect ^ \^f\ X 



42 cm-W"^ and |£p| = 0.0513 



eV Vsd 



+ 0.2676 eV for device A, 



then we can have a corrected G mode energy, Qg_t = ~ 
{AQG_in + AQG_defect)y ^^d we plot the new values in Fig. 4a as blue 
circle points. The new, after correction, G mode energies are fitted 
very well with a 2""^ order polynomial function. 



^G_r = cii V\ + a2 Vsd + , 



(2) 



where a^, a2 and are fitting parameters. For device A, we get = 
-0.0474 cm-^V-^ a2 = -1.6901 cm-'Y-\ and = 1580 cm''. 
We have not found any report about the coefficients of the graphene 
G mode energy as a function of the temperature for a temperature 
range up to 2500 K. H. Herchen and M. Cappelh reported that a first 
order Raman mode energy of diamond has a 2""^ order polynomial 
dependence on temperature up to 2000 K^^. So, it is reasonable to 
assume, for G mode in graphene, that 



(3) 



where hi, h2 and are graphene dependent only constants. The 
temperature of the device appears a linear relationship with Vsd, 
and we define 



T = ^Vsd + Tj,T. 



(4) 



where ^5 is a device and location dependent constant and Tj^r is the 
room temperature, 298.15 K. Combining equations (2), (3) and (4), 
we will have 



ai = bil]\ a2 = 2pTRTbi + b2P. 



(5) 



I. Calizo et al. reported that the graphene G mode energy linear 
coefficient with temperature is —0.016 cm"VK for a temperature 
range 83 K-373 K, and Qg_t = 1584 cm"^ for T = 0 Because 
bi,b2 and bs are graphene dependent only constants, we have 

AQg_t 



AT 



173 = 1584 cm 



-^b2= -0.016 cm~^/K, 



(6) 



in which 228 K is the average of 83 K and 373 K. Combining equa- 
tions (5) and (6) with knowing values for and a2 from fitting of 
device A, we can solve bi = -4.596 X 10"^ cm~^K~^, ^2 = 
-0.0139 cm-'K-\ b^ = 1584 cm-\ and P = 101.5 Wfor the mea- 
surement point in device A. Finally, we combine equation (1) and (4), 
and have \Ep\ =yT to describe a graphene device intrinsic Fermi level 
as a function of the temperature, where y = a/ip. From device A, our 
measurement gives out a graphene device intrinsic Fermi level, \Ep\ 
= 2.53 X 10"^ eVR-'-Tor \Ep\ = 2.93 ksT. Using graphene Hvp = 
1.16 X 10"^^ J-m, we have a graphene device intrinsic carrier density, 
Hir, = 3.87 X 10' cm-'K-'-r. 

In Fig. 4b, the square points are the measured G mode positions of 
device B as a function of Vsd- The circle points are the corrected G 
mode energies using the same correction method for device A, but 
.eV Vsd 



with \Ef\= 0.0565 



V 



+ 0.2676 eV. For the p value, it is invalid 



to process the limited data from device B with the same method for 
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Figure 5 | Device I-V curves, temperatures, and carrier drift velocities, (a), Device A, the defect level is 4 X 10^^ cm"^. (b), Device B. Red solid 
square points and black solid round points are the I-V data during Vsd downward and upward scan, respectively, using left y axis. The current values are 
the averages of the currents before and after the Raman measurements, the differences of the two are plotted as error bars. Red hollow diamond 
points and black hollow star points are calculated carrier drift velocities during Vsd downward and upward scan, respectively, using right y axis. Only the 
data points with carrier density /7 > 1.4 X 10^^ m"^ during the upward scan are plotted in order to control the error in drift velocities due to the 
uncertainty in Eg. The top axis plots the device temperature, which is calculated from Vsd using equation (4). 



device A. Luckily, ai is a fitting parameter independent from the 
values of a and Eg. In order to prove this point, we did the 2nd order 
polynomial fittings for both measured G mode positions and cor- 
rected G mode energies, and plotted the fitting curves as a dash line 
and a solid line, respectively, in Fig 4b. Both fittings give Ui = 
— 0.0603 cm"^V~^. Because bi is a graphene dependent only con- 
stant and should be same for all graphene, we can use the = 
—4.596 X 10"^ cm~^K~^ value from device A, for device B. Using 
equation (5), ai = hiP, we can get = 114.5 K/V for device B. 
Together with a ~ 0.0566 eV/V for device B, we get y = 2.47 X 
10"^ eV/K from device B. This y value of device B is consistent with 
the y value from device A. Such an agreement between device A and 
device B supports the validity of our method. 

Graphene shows a great potential as a material for high temper- 
ature devices, compared to traditional semiconductors. High tem- 
perature devices are demanded for military, astronautic, and heavy 
industry applications, where high power density is required or high 
temperature environment can't be avoided. In order to make a non- 
linear semiconductor device (p-n junction or transistor) functioning, 
the chemically doped carriers or the gating generated carriers have to 
be the controlling factor within the semiconductor material. 
However, when the device temperature increases, the semi- 
conductor's intrinsic carrier density will increase. Once the temper- 
ature reaches a certain point that, the semiconductor's intrinsic 
carrier surpasses the chemical doping or gating effect and dominates 
the device's response, the nonlinear device will simply fail its func- 
tion. For instance, a Si device's intrinsic carrier density will double if 
temperature increases 8 K from room temperature, the maximum 
operation temperature for a Si device is 250°C; a Ge device's intrinsic 
carrier density will double if temperature increases 12 K from room 
temperature, the maximum operation temperature for a Ge device is 
100 °C^°. Normally for a semiconductor, a larger band gap gives a 
higher operation temperature; however, a larger band gap also 
increases the difficulty of carrier density tuning. Graphene is semi- 
metallic, not a semiconductor. By a doping or field effect, graphene 
can be made into a p-n junction or a transistor as well. Running our 
numbers for graphene, we find that it will take 125 K for graphene 
intrinsic carriers to double from room temperature. At 1500 K, gra- 
phene intrinsic carrier density is 8.7 X 10^^ cm"^. Multiple research 
reports ^^'^^'^^ have demonstrated that a top gating method or a chem- 
ical doping can tune the graphene carrier density to a level above ±3 
X 10^^ cm"^. Therefore, we foresee that a well-fabricated top gated 



graphene transistor or a chemical doped graphene p-n junction will 
be still well functioning at a temperature above 1500 K. A lot of 
things melt at 1500 K, for example gold melting temperature is about 
1340 K. So, such a graphene transistor or graphene p-n junction will 
never fail due to rising intrinsic carrier density before other parts start 
to melt. Above analyses reveal the great potential of graphene as the 
material for high temperature devices and point out a new research 
direction for graphene applications. 

Finally, we will show and discuss carrier drift velocity in graphene 
up to 2400 K. Fig. 5 shows the device I-V curves, temperatures and 
carrier drift velocities of device A and device B. Carrier drift velocity, 
V, equates to I/(qWn), where I is the current, q is the electron charge, 
W is the device width, n is the carrier density. The device B broke 
down after Vsd > 17 V and the device A was also on the brink of 
breakdown Vsd = 21 V and showed clear runaway effect, in which 
the current keeps increasing and feeds itself with a constant driving 
Vsd- The maximum temperature is about 2400 K for device A and 
2200 K for device B, those temperatures are consistent with the 
theoretically calculated graphene breakdown temperature, 2230 ± 
700 K, in vaccum^. The drift velocities from device A, Vsd downward 
have the best accuracy compared to the values from upward scans, 
and map out a very nice three stage behavior under high field. In the 
first stage, the drift velocity increases linearly with the field. Ohm's 
law holds perfectly. Then with a further increased field, lattice scat- 
tering dominates, the drift velocity still increases but with a reducing 
rate. Ohm's law still works but with a reducing conductivity. In the 
last stage, after the field reaches a certain level, intrinsic carrier dom- 
inates, the drift velocity stops to increase. Traditionally, the drift 
velocity at this point is called the saturation drift velocity, and this 
saturation drift velocity stays the same for traditional semi- 
conductors if the field increases further. For example, for Si, it is 
about 1 X 10^ m/s. Our graphene sample shows a linear drop of 
"saturation" drift velocity with the increasing temperature. This 
reducing saturation drift velocity of graphene has been predicted 
by theoretical calculation considering the degeneracy of carriers 
and the Pauli exclusion principle^^. The same calculation gives a 
saturation drift velocity value, 2 X 10^ m/s, for an ideal graphene 
with a doping level 10^^ cm"^ at 300 K. The maximum drift velocity 
of our graphene with a defect level of 4 X 10^^ cm"^ is about 4 X 
10^ m/s at 1000 K. Our value is not far from the calculation value 
after considering the temperature difference and the scattering from 
the defects. The declining saturation drift velocity in graphene helps 
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delay the current runaway point, where the current keeps increasing 
on its own with a constant driving voltage, to an even higher tem- 
perature; but it also limits the ability of graphene in delivering large 
current. Another thing we want to point out is that Ohm's law fails 
totally in the intrinsic carrier dominated third stage. The electrons 
and the lattice are in a non-equilibrium state, in which the energy 
electrons taking from the field is not fully dumped to the lattice. This 
is why the device temperature has a linear relationship with V, not 
the IV. The difference in drift velocity for a same Vsd between down- 
ward and upward scan is due to different effective electron tempera- 
tures, despite a similar lattice temperature. During downward scan a 
graphene device appears a smaller contact resistant (or barrier) than 
the upward scan^^, therefore, downward scan data has a larger cur- 
rent, electrons carrying more energy and moving faster, or a higher 
effective electron temperature than the upward scan data. 

In summary, we use Raman method to quantitatively measure the 
graphene intrinsic carrier density, drift velocity and G mode phonon 
energy as functions of temperature up to 2500 K. We find that the 
intrinsic carrier density of graphene is an order of magnitude less 
sensitive to temperature than that of traditional semiconductor 
materials. We foresee graphene becomes the material for high tem- 
perature devices and a properly fabricated graphene transistor or p-n 
junction will not fail at all before other parts melt. We map out the 
three stage behavior of carrier drift velocity as a function of the field 
or temperature, and graphene shows a linear decline in "saturation" 
drift velocity in contrast to a constant saturation drift velocity in 
traditional semiconductor materials. We also find the G mode pho- 
non intrinsic energy coefficients with temperature, — ~ 4.596 X 
10"" cm-'R-^T + -0.0139 cm-'R-'T + 1584 cm-\ All of above 
reveal the great potential of graphene as a foundational material for 
high temperature devices, identify a new research direction for gra- 
phene applications, and are vital in understanding the physics of 
graphene under high field or high temperature. 
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